From non equilibrium quantum Brownian motion to impurity dynamics in 
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Impurity motion in one dimensional ultra cold quantum liquids confined in an optical trap has 
attracted much interest recently. As a step towards its full understanding, we construct a generating 
functional from which we derive the position non equilibrium correlation function of a quantum 
Brownian particle with general Gaussian non-factorizing initial conditions. We investigate the slow 
dynamics of a particle confined in a harmonic potential after a position measurement; the rapid 
relaxation of a particle trapped in a harmonic potential after a quantum quench realized as a 
sudden change in the potential parameters; and the evolution of an impurity in contact with a one 
dimensional bosonic quantum gas. We argue that such an impurity-Luttinger liquid system, which 
has been recently realized experimentally, admits a simple modeling as quantum Brownian motion 
in a super Ohmic bath. 



C/3 

t>JO. 



o 

> 

m 

in 

rn 

O 



X 



I. INTRODUCTION 



Quantum Brownian motion has been the starting point 
for the understanding of more complex dissipative quan- 
tum systems [l[ . Applications to quantum tunnel junc- 
tions y, dissipative two-state systems [3| and reaction- 
rate theory [J] are just a few among many. In its simplest 
form, as proposed in the founding papers 0-0] , the envi- 
ronment induced dissipation is modeled by an ensemble 
of quantum harmonic oscillators linearly coupled to the 
particle of interest. So far, in most studies of the dis- 
sipativ e dy namics of a harmonically confined [y, |8| or a 
free [y, li3l quantum particle, the quantity of interest has 
been the reduced density matrix that is obtained by trac- 
ing away the bath degrees of freedom in the density ma- 
trix of the coupled system. For generic initial conditions 
this quantity has been obtained with the help of func- 
tional integral methods @, [Til. An alternative simpler, 
though in general only approximate, description of the re- 
duced density matrix is given by a master equation. For 
factorizing initial conditions [l2| and thermalized initial 
conditions [13| an exact master equation can be obtained. 
However, it is also known that there cannot be a master 
equation - in the form of a partial differential eq uation 
local in time - for arbitrary initial conditions |13| . The 
alternative quantum Langeyin approach |14l | extensively 
used in quantum optics [15|] is not sufficiently powerful 
either, for quite the same reason: only a few special ini- 
tial conditions can be successfully treated within this ap- 
proach and the quantum noise statistics are not tractable 
in the generic case. 



Non equilibrium correlation functions have been re- 
cently observed in cold atoms experiments on the dy- 
namics of an impurity atom moving in a one dimensional 
(ID) quantum liquid J^ J8j. Both the impurity and 
the quantum liquid are confined in an optical harmonic 
trap so that the impurity motion resembles the dynam- 
ics of a damped quantum harmonic oscillator. In [l9[ 
the authors attempted to describe the impurity dynam- 
ics within the Gross-Pitaevskii approach at zero temper- 



ature. In this paper we will follow an alternative way by 
applying quantum Brownian motion theory to the im- 
purity problem. More specifically, we will generalize the 
path integral formalism in [0] to the language of gener- 
ating junctionals commonly used in quantum field the- 
ory. From the generating functional we will easily deduce 
the non- equilibrium correlation Junctions for generic non- 
factorizing Gaussian initial conditions (for a stochastic 
description of open quantum systems via a generating 
functional, see (20[) that cannot be obtained within the 
density matrix formalism. We will use this technique to 
treat three problems: the relaxation dynamics of a par- 
ticle confined in a harmonic potential after a position 
measurement performed at the "initial time" ; the relax- 
ation dynamics of a particle after an abrupt change in the 
parameters of the confining potential; the motion of an 
impurity [iSl [2l| - |24| in contact with a ID quantum gas all 
confined in a harmonic potential. In the latter case, we 
will explicitly compare our theoretical findings to recent 
experimental results |16| on ultra cold quantum gases. 
We will argue that, although the impurity-Luttinger liq- 
uid system is described by the Frohlich polaron Hamil- 
tonian, many aspects of the impurity dynamics can be 
understood in the framework of quantum Brownian mo- 
tion. 

More precisely, the paper is organized as follows. In 
Sec. |II] we present the model and we review the main 
results obtained in Q. In Sec. [ITT] we employ path inte- 
gral methods to derive the generating functional of out 
of equilibrium correlations. Our results cover both fac- 
torized and non-factorizcd Gaussian initial conditions as 
well as the effects of an initial position measurement per- 
formed on the particle. In Sec.|IV]we study the equilibra- 
tion processes after an initial position measurement and 
after a quench in the harmonic potential and we derive 
the equilibration times for low and high bath temper- 
atures. In Sec. |V] we apply our formalism to impurity 
motion in a ID quantum gas described by the Luttinger 
theory. To keep the discussion simple, we choose to use a 
simplified modeling of the experiment in which we neglect 



polaronic effects [25|, |26| as well as the possible renormal- 
ization of the external potential [1^1 ■ These subtle ef- 
fects will be analyzed elsewhere. The Luttinger liquid is 
found to behave as an exotic quantum bath of harmonic 
oscillators with a highly non Ohmic spectral density and 
non-linearly coupled to the particle. This is shown to 
lead to the curious behavior that the oscillator frequency 
can increase upon increasing the coupling constant be- 
tween the "bath" and the impurity in strong contrast to 
the behavior of an Ohmic damped oscillator. We fur- 
ther calculate the non equilibrium equal time correlation 
function (the variance of the position) and we compare 
our theoretical results to experimental evidence. In the 
last section, Sec. IVH we conclude and we present further 
possible applications of our work. 



II. THE MODEL 

We study the evolution of a particle of mass M evolving 
in a (possibly time-dependent) potential V{q;t) where q 
is the position operator. The Brownian motion stems 
from its interaction with a quantum heat bath which is 
usually modeled by an infinite set of harmonic oscilla- 
tors linearly coupled to the position operator q. The 
full system is then described by the Hamiltonian T-L = 
lis + Hb + TisB, with 



all possible observablcs A is governed by 
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(5) 



with T the time-ordering operator. By introducing the 
density matrix of the initial state po the A^-time average 
of a set {Ai} oi N operators, is 



{ANitN)AN-litN-l) ■ ■ ■ Ai{ti)) 



Tr 



ANitN) ■ ■ ■ Ai{ti)po 



(6) 



where we took the product of the AiS to be time ordered 
(with t]s[ > ijv-i > . . . > ii) so that we can more easily 
make the connection between Eq. ^ and its path inte- 
gral representation. We assumed that Trpo is normalized 
to one. Note that for a generic initial matrix po this, as 
well as any other, correlation function is not necessarily 
stationary, i.e., it may depend on the N times explicitly. 
In all cases the model has to be supplemented by in- 
formation on the initial condition of the coupled system. 
These are incorporated in the initial density matrix po- 
Equilibrium dynamics can be studied by choosing po to 
be the Boltzmann weight, that is 



ns[q,p]^j^ + V{q;t)-H{t)q, 



TiBiiXn, Pn}] = X! 



2mr, 






TisB[{Xn},q] = -q^CnXn + f^ 1^:;;^^, 
n=l 
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(1) 
(2) 

(3) 



p is the momentum operator of the particle. The last 
term in Eq. ([l} introduces a time-dependent source H{t), 
a c-number, that couples linearly to the particle's posi- 
tion q. Xn and p„ are the position and momentum op- 
erators of the n-th harmonic oscillator, with mass and 
frequency 7ti„ and w„, respectively. c„ is the coupling 
strength between the particle and the n-th oscillator's 
position. The last term in Eq. ([3]) compensates for the 
bath-induced renormalization of the potential. Indeed, 
the sum of Eqs. © and ^ can be rewritten as 



Hi 



nsB = Y: 



PJ 
2mr, 



mnUjl 



m„w2 



9? , (4) 



which shows the absence of any drift force induced by the 
bath and ensures that V{q, t) corresponds to the physical 
potential right from the start. The model Hamiltonian 
Eqs. ©-(H]) has been widely used in the literature as a 
generic model for the dissipative dynamics of a quantum 
particle P,0,il- 

In the Heisenberg representation the time evolution of 



Po = exp(-/3-H) 



(7) 



where % is the full coupled Hamiltonian and the nor- 
malization constant has been ignored. This truly equi- 
librium density matrix has to be distinguished from 
Po = e~^^^ ® q-0'Hb ^ g^ (jg^gg jjj which each component 
of the "universe" (the whole particle-bath system) is in 
equilibrium on its own at the same temperature. This 
subtle point is often overlooked in the literature. 

Non equilibrium dynamics can be studied whenever the 
initial density matrix is not of the form in Eq. ([7]) . The 
simplest choice is an initial product state for which the 
initial density matrix factorizcs into two contributions 
pso and pBo which solely depend on particle and bath 
variables, respectively: 



Po = Pso (X) PBO 



(8) 



Brownian motion [l|, |27| as well as the dynamics of more 
complex macroscopic systems [2a - [3a | with a factorized 
initial density matrix have been studied in a variety of 
physical situations. However, in many cases it is not 
appropriate to assume Eq. ([8]) since one has no command 
over the bath and it is impossible to "switch it on and 
off" at will. In addition, with recent developments in 
cold atom experiments, new classes of initial conditions 
become of relevance. 

The first one covers all situations in which the particle 
is in equilibrium in a potential and either it is released or 
the potential is suddenly modified at t = 0. In this case 



Eq. ([7]) holds with T-L replaced by Ho (7^ "H) describing 
the initial state. 

The second class concerns all situations in which the 
position of the free particle is measured at f = 0. This 
procedure projects the initial density matrix onto the 
quantum states of the measurement outcome. We fo- 
cus on the case where no quantum quench is performed 
in addition to the position measurement so that "Hq = ^■ 
If the position is exactly determined at t = the initial 
density matrix is 

Po = tl{q)e-"^il{q) with Ii{q) = \q){q\ (9) 

the projection operator onto the state \q). If, instead, we 
take the measured position of the particle to be Gaussian 
distributed around g„i the projection operator takes the 
form 



n(g„^) 



dg e 



l'Z)('Zl, 



(10) 



where A measures the uncertainty of the particle's po- 
sition at t = 0. Once again we neglected the irrelevant 
normalization factor. 

A third important class of initial conditions are the 
factorized density operators, see Eq. (|8]), in which the 
initial state of the system is a pure state. Since any 
state can be expanded in terms of displaced Gaussians 
(or coherent states) it suffices to consider initial states of 
the form 



Pso = \'>P){v\ 



where 



^{q) 



4a2 



and ^{q') ~ e 



{<,'- 



(11) 



(12) 



to cover the whole class of initially factorized pure states. 
In this article we derive a generating functional that al- 
lows us to obtain the A^-time correlators for these types 
of initial conditions. We are mainly interested in the evo- 
lution and averages of the particle's position observables 
for which A = A[q] with some function A depending on 
the position q of the particle. Note that due to the cou- 
pling of the bath to the particle's position the momentum 
dynamics follow from the Heisenberg equation Mdtq = p. 
Therefore, by focusing on the particle position operator 
we simultaneously describe the dynamics of the parti- 
cle's momentum. While in Q the authors derived an 
explicit expression for the equilibrium correlation func- 
tions the generating functional will allow us to go beyond 
the equilibrium case. 



III. THE GENERATING FUNCTIONAL 

In this section we derive the generating functional of 
all (non-equilibrium) correlation functions. This goes 
beyond the analysis in [9| where explicit expressions 
for equilibrium correlation functions were given. More 
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precisely, we derive a functional of two time-dependent 
sources J[F, G] such that the two-time correlation is 
given by 



(13) 

F,G=0 

and similarly for higher order correlations. 

We obtain the path integral formulation of the gener- 
ating functional by making use of the coherent states of 
the bath variables |{Cn,/}) which are defined in App. \K\ 
The ensuing functional integration includes paths over 
particle and bath variables. Since we are not interested 
in the degrees of freedom of the bath, we average over 
all bath variables to find a "reduced action" that only 
depends on the particle variables. In the special cases 
discussed below (e.g., harmonic potential) the remaining 
path-integrals can also be performed and the functional 
J can be fully determined. In this section we sketch 
all steps in the derivation. Further technical details are 
reported in App. [X] The reader who is not interested 
in these technical details can jump directly to Eq. ([M)) 
where its rather lengthy final expression is given. 



A. The density matrix 

In terms of the product states \q, {inj}) between the 
particle and the bath eigenstates the matrix elements of 
the time evolution operator read 

{<lfAUf}\fe-i^o'^"^^'\.AinA) (14) 
and of its Hermitian conjugate 

/c*(q;,{e;/};<zUc;j;0^ 

T^ is the anti-chronological time ordering operator. The 
elements of the time-dependent density matrix, p{t) = 
/Cpo/C*, are given by 

w(<z/,{en,/};g/,{C,/};i) ^ (9/,{en,/}| m WfACj}) 

= / dgrjdg-d^id^- IC{qf , {^nj}] Qi, {^n,iht) 

X W(g„{C„,J;<z,',{C,:}) IC^q'fACjhQlienA-.t) 

(16) 

where the matrix elements of the initial density matrix 
have been denoted by 

W{q^,{^n,^hq^AC.^}) = {q.,{Cn4\ PO H, {C.^}) (17) 

and we used the short-hand notation 



d6 = n e-«".-«-de:,,dc„ 



(18) 



The path integral representations of K, and /C* are 



with 



where we made clear with the superscripts ^ and ^ which 
paths belong to /C and /C*, respectively. The functional 
integration measures are defined in App. |^ 



B. Reduced density matrix for a system initially 
coupled to an equilibrium bath 



The time-dependent density matrix in Eq. (J16p still 
contains information about the degrees of freedom of the 
bath which we are not interested in. Therefore, we aver- 
age (trace) over all the bath variables to find a reduced 
density matrix that depends only on the particle vari- 
ables and the external sources. 

We are interested in a system that is initially coupled 
to an equilibrium bath. Therefore, 



'H-Bo — T~Li 



n 



SBO 



^n 



SB 



(19) 



where all initial Hamiltonians are labeled with a subscript 
Q. At this point it is not necessary to make Hso explicit 
since this term involves only particle variables that are 
not affected by the trace over the bath variables. The 
matrix element of the initial density operator Eq. (jl7|) in 
Eq. (J16p can be represented by an imaginary time path 
integral 

where the initial action So is in general different from S 
reflecting the fact that T-Lq ^ H. The reduced density 
matrix can now be recast as 

Wiqf, q'f\t) = j dqMvq"T^q^T^q~ 

where J-'[(7+, g~,(7°] is the "influence functional" that de- 
pends only on the particle variables, as also do Ss and 
Sso- The path integral runs over all paths with 

q+{t) = qf, g+(0)=(7, , q- {t) ^ q'j , 
q-{0)=ql, g°(/3?i) = g-(0), gO(0) = g+(0), 

which is the reason for the name "closed-time path inte- 
gral". It is convenient to introduce the linear combina- 
tions 



X=iq+ + q-)/2 



and 



X = q 



q 



(20) 



The calculation of J^ can be found in App. lAll or in Q; 
the result reads 



F[XtX, q^] = e 



ir-i>[x,x,q°\ 



(21) 



$[a:,a;,q"] 



ph 



drdcr k{T-a)q^{T)q^{a) 



+ At ds K*{s~iT)q°{T)x{s) 

Jo JQ 

if* 
— — I dsdu K ji{s — u)x{s)x{u) 

2 Jo 

-M / dsi(s)— / Au -i{s ~ u)x(u) . (22) 
Jo ds 7o 

The kernels K{d), 7(i) and k{t) are defined in 
Eqs. (|A3ip . (|A42p and (jA43|) . respectively. Kfi denotes 
the real part of K. Note that $[x, x, g"] depends on the 
fixed "end-points" t and Ph of the closed-time path. 

Expected values evaluated at different times are now 
expressed in terms of a path-integral over g°, x and x 
with an effective action S, 



dxfdxidxi / 'Dx'DxDq 



^^^[x,x,q°] 



(23) 



where T,[x,x,q^] is given by 
^x, X, q"] + iSso[q°] + Ss[x + x/2] - Ss[x - x/2] 



^x,x,q°] +i 
t 



ds 



Mxx 



I3h 



Vh 



dr 



^iqr + Vo{q') 



X 

x+-;s 



Vh' 



x 



(24) 

We introduced the initial mass Mq of the particle and the 
initial potential Vq that are in general different from the 
"bulk" mass AI and potential V. This allows for quenches 
in these parameters. Note that the case in which the 
initial state is a pure state [e.g. Eq. PTj) ] can be easily 
recovered by setting Mq = and Vg = or, equivalently, 
by noting that the path 5° shrinks identically to zero 
(since there is no initial Hamiltonian for this simple type 
of initial condition). 

The superscripts in the path integral in Eq. (j23p re- 
mind us of the constraint that the paths arc subject to. 
One has 



x(0) 

Xit): 



Xi, x{0) 
x^ , 
,0 



f, ^0(0)=g+(0)=x, + ^, 



(25) 



x(i)=0, q'>il3h)=q-{0)^x. 



~2 



Note that due to the periodic boundary conditions of the 
trace Xf ~ q'^{t) — q~{t) = 0. 

C. Generic Gaussian initial conditions 

It is very easy to include the change of $[x,x,g°] in- 
duced by the initial position measurement in Eq. ^ . By 



using the explicit Gaussian form of the projector n(gm) 
[see Eq. PH)) ] the dependence on the initial measurement 
can be simply incorporated in ^[x, x, g°] by an additional 
term of the form 



sources reads 



4A2 


(X, + - - 


ih f 




2A2 



<7m)^ + iXi- -^ - QmY 



{Xi ~ g„j)^ + 



(26) 



In the limit of strong uncertainty A — > oo the effect of 
the initial measurement is blurred. 

In order to recover the case where the initial state of 
the system is pure and decouples from the environment 
[which corresponds to the factorized initial density ma- 
trix with psQ given by Eq. (|lip ] the action in Eq. ([M]) 
has to be supplemented by 



ih 
4A2 



[Xi 



ih 
2A2 



with the notation 



y - Qrnf + {X^ - 


X^ 
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-2 -2 

, 2^J , 2 , ^7n 


-2x 



q'mf 



i-Xi (.. ^n 



{q-m + (im)l'2 and Xm = (g« 



(27) 



(28) 



Since Eq. ([26|) is a special case of Eq. (f27| corresponding 
to qm = q'm (or Xm = q-m and Xm = 0) we will work with 
the latter in the following. The relevant cases can then 
be selected by taking simple limits. 

In the following expressions we will write only the 
terms that depend on Xi or Xi since the ones depending 
on Xm and Xm contribute only to an overall constant. 



D. The sources 

The source term appears as J dt'H{t')q^{t') in K, [see 
Eq. (mi)] and as - J dt'H{t')q-{t') in /C* [see Eq. p^])]. 
For convenience, we distinguished the function existing 
on the positive running branch of the closed time contour, 
which we still call H{t), from the one existing on the 
negative running branch of the same contour, which we 
call H'{t). This implies that the potentials in Eq. ([24| arc 
given by Vw (y) ^V{y)~ Hy and Vh' [y] ^V{y)- H'y. 

After the transformation of variables in Eq. ([20]) we ob- 
tain two external time-dependent sources F{s) = [H{s) + 
H'{s)]/2 and G{s) = [H{s)-H'{s)] which couple linearly 
to the variables x(s) and a;(s), respectively. All correla- 
tion functions can be computed from the generating func- 
tional J^[F, G] as derivatives of J' with respect to i^ or G 
evaluated at i^ = G = 0. A physical force is represented 
by H{s) = H'{s), that is by F{s) ^ and G{s) = 0. 
Therefore, the linear response of the mean value Eq. ([25)1 
to an external force can be obtained for F{s) ^ 0. 

The generating functional, that is to say, the trace over 
the reduced density matrix in the presence of the external 



gJ[F,G] r^ J dxidx.dxf J' T>xT>xT>q° ei^^"^'^' 



q ,Xi,Xf,Xi 



where the path integral is subject to the constraints in 
Eqs. ([25]) . The overall normalization factor depends on t, 
(3 and all parameters in the model but not on the fields. 
We can now write 
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(29) 



and 



mm')) = 

h 5 h 
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.J[F,G] 



F,G=0 



(30) 



and all other correlation functions can be obtained in 
a similar way by noting that q^{t) = qf = Xf and 
q^{t') = x{t) + x{t)/2. At this point it has become ob- 
vious why two sources are needed in order to obtain all 
non-equilibrium correlation functions. 



E. The harmonic case 

To go further we restrict ourselves to the study of a 
quantum Brownian particle in a harmonic potential for 
which 



- V{x + x/2; s) + V{x - x/2; s) = --Mw'^xx (31) 



and 



^o(g") 



-A/o^o(9")' 



(32) 



The choice of a quadratic potential renders the problem 
analytically solvable. The generating functional can be 
calculated by simply evaluating the action on its minimal 
action path (over the initial condition branch and the 
time-dependent branches) as Gaussian fluctuations yield 
only pre-factors that are independent of the sources and 
can be determined at the end of the calculation from the 
normalization of the density matrix. Note that, although 
both initial and bulk potentials are harmonic, they are 
not necessarily the same thus allowing for the study of 
quantum quenches. 



F. Integration over the initial condition 

We first treat the contribution of the initial condition 
path g" in Eq. (p4)) . The equation of motion for (f can 
be easily obtained from Eq. (P^ : 

Mofir) - / da k{T - a)q°{a) - A/ow^^^r) 
Jo 

= -i ds K*{s - iT)x{s) , (33) 

Jo 



with the fixed end-points q^{0) = qf and q^{l3h) = q~ . 
As the q° path is part of the whole closed-time path it 
implicitly depends on the fixed end-time t as well. In Q 
one can find a detailed analysis of this equation of motion 



which uses a Fourier expansion of the path q^ (r) on the 
interval [0,(3h]. By using the results found therein we 
obtain 



Zj^X^ X^ Xj, X j^ ^ Xi\ 
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+ T7Fis)^is) + T7Gis)x{s)-x{s)— / dwy^s - u)x{u) 



M 



M 



ds 



(34) 



We introduced the complex "force" 

F{s) ^ F{s) + x.Ciis) ~ ix,C2is) 
with the functions Ci and C2 



(35) 
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ptiK 
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^ J2 UkVkfk{s) 



C2{S) 

Dfl. 
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The constants A and Q are given by 



(36) 



A:— — 00 k— — oo 

(37) 

withufc = {LO^ + iyl + Ck)-^ Vk = 27rfc/^?i, Cfe = [Af7(0)- 
gfc(0)]/Afo [for the definition of 7(i) see Eq. ^ below] 
and 

/>00 1 rj 

9k{s)^ — S'(w)^- — 2Cos(a;s), (38) 



fkis) 



duj 

— S'(w, 



2i^fe 



■ sin(a;s) 



(39) 



/O " "^ ' ''fc 

where S{uj) is the spectral density of the bath. The two- 
time function R' reads 

R'{s, u) = R{s, u) + MoKr{s - u) , 
i?(s,u)==-ACi(s)Ci(u) 

^ Ufc[5fe(s).9;.(ii)-A(s)A(ii)] , (40) 



with 



im 



Kr{s -u)^ —Y^ gk{s - u) 



(41) 



the real part of the kernel K. The time-dependent bath 
kernel 7(5) is given by [see Eq. (|A42p ] 

duj S{uj) 



M Jo TT 



■ cos(ws) 



(42) 



The functions Ci and C2 as well as the kernel R{s, u) are 
not to be confused with the correlation functions and the 
linear response function that will be denoted by C and 
TZ, respectively. 



G. Real time minimal action paths with external 
sources 



The equations of motion for x{s) and x(s) read 



i{s) 



d 

ds JO 

F{s) 



du 7(5 — u)x{u) + uj x{s) 



du R'{s,u)x{u) , 



^(^) - d^ 



M MMa Jo 

du 7(u — s)x{u) + ur'x{s) = 



M 



(43) 
(44) 



The action E evaluated along the minimal action paths 
can be determined by inserting the solutions to Eqs. ((32]) 
and (j44p into Eq. (p4)) . However, the authors of Q noted 
a simplification of the calculation which can be general- 
ized to our case where the source G{s) is also present (in 
[9| no external source for x was used). The idea is the 
following. After a partial integration in the second line 
of Eq. p4p the action E takes the form 



dsdw R'{s,u)x{s)x{u) 



Zjj^X, X, Xi^ ^fi ^i 
— MXiXi 



2M. 



Jo 



+ / ds G'(s)a;(s) + border(a;i,a;/,a:i) (45) 

Jo 



when evaluated along the minimal action paths deter- 
mined by Eqs. (|43|) and (J44|) . where we used the bound- 
ary condition Xf — 0. Here, border(a;i,a;/, Si) stands for 
all border terms in Eq. ([M]). 

On the other hand, one can split the force Eq. ((35)) 
into its real and imaginary parts F{s) = Fii{s) + iFi(s). 
Then, the minimal action path x{s) splits into x{s) = 
a^fl(s) + 'ixj{s), where xj{s) satisfies the boundary con- 
ditions xj{Q) = xi{t) = 0. The trick is to show now 
that one can simply focus on the real part xr{s) of the 
minimal action path in order to obtain the complete sta- 
tionary phase action. Indeed, if we evaluate the action 



Eq. ((34)) only along the minimal xji{s) and x{s) we obtain and 



T,[xR,x,Xi,Xf,Xi] = -MxtXii^i+ / ds G{s)xr{s) 



+ ds x{s) 



Fiis) + 



2A/( 



Jo 



du R (s, u)x(u) 



-hoTdeT{xi, X f , Xi) 



(46) 



X(s) — — -^ Xi - 

M g+{t) 



M 



du g+{u — s)G(u) 



du g+{u)G{u) , (48) 



where we used the fact that xr{s) satisfies the real part 
of Eq. (1321) • We now want to show that Eqs. (gB]) and (|i5)) 

are indeed equal. With the help of the imaginary part of where g+{t) is a propagator that in Laplace-transform 

reads 



Eq. (|43)) and the equation of motion (|44| we can easily 
prove by integration by parts that 



dsx{s) Fi(s) + —— / du R'{s,u)x{u) 
Mq Jo 

= —MxiXji + ds G{s)xi{s) , 



and by using this identity in Eq. (j46|) we recover Eq. (|45p . 
Therefore, the right-hand-side (rhs) of Eq. (|45)) and the 
rhs of Eq. (^5)) coincide. It is sufficient to evaluate the 
action Eq. (|34p along the real component xii{s) that sat- 
isfies a much simpler equation than x{s). 

In terms of the end points Xi, Xf and Xi, the solutions 
to the real parts of Eqs. (|43]) and (04]) read 

Xr{s) = ^ ^ Xf 



1 



O.M-'^M.[ 



M 



— / du g+{s - u)Fr{u) 



1 G+js) 
MG+{t)Jo 



G+iX) 



A2 + A7(A)+w2 



(49) 



From Eq. (j47|) we immediately find, by using the bound- 
ary conditions ^+(0) = 0, ^+(0) = 1 and ^+(0) = 0, 



XR,; 



;^f 



G+{t) 



G+(t) ' G+{t) 



Xi 



du Q+{t - u)Fr{u) 



(50) 



du G+{t — u)Fr{u) (47) external sources F and G: 



1 1 

G^MJo 



Inserting the solutions to Eqs. (J47l) and (|48l), and Eq. ([50]) 
into Eq. (|46p we find an effective action I][xi, x/, Xi] that 
depends only on the end-points, x^, Xf and Xi, and the 



ZjyX^ , X f , Xi\ 



iMo 

.2 - 



iMn 



-XiXj^ 



MxiXf- 



G+{t) 



'G^i^)-MlGAt) 



G+it)- 



-Xi I ds G{s) 
'o 

-^ I* dsdu^±^^g+is)Gis)FRiu) 



.._ g+{t) X, 
'g+{t) G+{t)Jo 



du g+{t - u)Fr{u) 



' r^ c ( \ 

dsdu R\s,u)x{s)x{u)+Xf ds-t-— 



— / ds / du g+{s - u)G{s)Fr{u) 



Gis) 



G+{t) 



(51) 



with 



1 



1 1 



A' A 
and 



r2' = fi 



4e2 



^2_MoA2 



(52) 



For the sake of a clear presentation we have not replaced 
x{s) and Fr by their corresponding expressions in terms 
of the end-points yet. 



H. Integration over the end-points 

In order to find the final expression for J'[F, G] we still 
have to integrate over the end-points Xi, Xf and Xj. Since 
the exponent Y,[xi,Xf,Xi\ is linear in Xf the integration 
over this variable generates a (5-function of the form 



Xi - — / ds g+{s)G{s) 



(53) 



up to a factor not depending on the end-points which, jj J du Q^{u — s)G{u) [the first and the third terms 
in combination with the integration over Xi, enforces a of the rhs of Eq. (pS)) cancel]. Moreover, after the inte- 
substitution of x.i by jj J^ds g+{s)G{s) and x{s) by gration over Xi the fifth and the last terms of the rhs of 

Eq. ([5T|) cancel, too. Finally, the Gaussian integral over 

Xi yields 



J\F,G] 



A' 



2hMo vjo 
A' '•* 



dsGis)g+is)] + 



Mq 



Xm / ds Q+{s)G{s) 



hMMo 7o 



t rs O'Af 

dsG{s)g+{s)l ds' du'g+{s'-u')G,{u')G{s')~^ 



t ^2 

ds G{s)g+{s) 



'iu I '^^ f '^'^ ^+^' " u)G{s)F{u) + J- /■ ds g+{s)G{s) I ds' f du'g+iu' ~ s' 



)G2{s')G{u') 



ds G{s)g+{s) + — ds du g+is - u)Gis)Giiu) 
^" Jo Jo 






2hMmoJo "^"io '^'' ' '^'^ ' ^l"'^+("-*)^"(^'«')^+("'-^')^(")^(^') 



(54) 



Here we introduced the kernel 



I. The correlation function 



The two-time correlation function Eq. ([5n| has two 
contributions. 



R"{s,s')^A'Gi{s)Ci{s')+R'{s,s') , (55) 



= C{t,t') + iA{t-t') . 



(56) 



The first term, the average of the anti-commutator or 
where we used Eq. (|A36p . Equation ([5^ is the central symmetrized contribution, is real and the second one, the 
result of the first part of this paper. It allows us to derive average of the commutator or anti-symmetrized contribu- 
all non-equilibrium correlation functions in a systematic tion, is imaginary and proportional to the linear response 
way. Direct applications of this method will be presented function, TZ{t,t'), as shown by the Kubo formula. For 
in Sees. IIVI and IVl generic Gaussian initial conditions, Eq. pH)) . one finds 



C{t,t') 



A'h 



Mo 
hfl'Mo 



g+{t)g+{t') + ^ j du g+it' - u)c,iu) + ^g+{t') J du g+{t - u)Gi{u) 



^, g+{t)g+it')- — 



h 



G+it) I Q+{t'-u)C2{u) + g+{t') [ g+{t-u)G2{u) 



A" 

H — i-x 



4 "^m 



A/ 



S+it) + T7 / du g+{t - u)Ci{ii) 



Q+it') + 



1 

M 



du g+{t' ~ u)Gi{u) 



^° -,xlg^{t)g^{t')+ ^ 



IGe^AP 



t' 



APMo Jo Jo 

r 



du / du' g+{t - u)R"{u,u')g+{t' - u') . 



(57) 



To represent a pure initial condition that is initially de- and R"{u,u') = MK{u — u'). A non-factorized initial 
coupled from the bath, as in Eq. (|lip . we set M = Mo, state with ~ say - an initial position measurement of 
A' = e^ and il' = l/4e^ as well as Ci(s) = 0, 6*2(5) = "width" e is obtained with Xm = Qm and Xm = 0. In the 



classical case the non equilibrium correlator of a Brown- 
ian particle is given by Eq. (|B3|) . 



IV. 



NON EQUILIBRIUM DYNAMICS AFTER 
QUANTUM QUENCHES. 



Quenches from a high temperature initial state have 
been extensively studied in the literature. They corre- 
spond to the case of a factorizing density matrix as in 
Eq. (|S]). In this section we will study the non-equilibrium 
dynamics of a quantum Brownian particle after quenches 
from different initial states. As already mentioned in the 
Introduction, there are two experimental scenarios that 
are of interest to us. In the first one the initial position 
of a Brownian particle trapped by a harmonic potential 
is measured at t = 0. We will study this case in the first 
part of this section by focusing on the two-time correla- 
tion function. The second scenario consists in a quench 
of the trapping potential, which in the case of a harmonic 
potential corresponds to an abrupt change in the trap- 

I 



ping frequency. Such a quench will be studied in the 
second part of this section. In both cases, we will derive 
the asymptotic equilibration behavior of the system in 
the presence of Ohnfic dissipation. 



A. A particle in a harmonic potential with an 
initial position measurement 

The general results in Sec. IIIII are here specialized to 
the case of a particle trapped in a harmonic potential, on 
which a position measurement is performed at f = 0. As 
we work with the same particle initially and subsequently, 
Mq = M, while loq = lj > 0. At f = a measurement of 
the particle position is performed with outcome q,n = 
and uncertainty A. The initial density matrix is given by 
Eq. ([9]) . Note that the particle is permanently coupled to 
the bath, hence the initial density matrix does not fac- 
torizc. Thus, our starting point is the general expression 
Eq. (|57p with M = Mq, lo ~ ujq and Xm ~ x„i = 0. The 
Laplace transform of the correlator reads 



c(A,«) - ^g+{x)g+{^) U'XK + ^xc,{n) + ^^^c,ix) 



n'Mg 



M 



C2{X) - —C2{n) + —.R"iX, k) ,(58) 



M 



A'P' 



r 



an expression that can be simplified by using the method 
explained in App. [Cl 

Introducing the function hk{X) — gk{X)/M 4- A, 
Eq. ([55)1 can be written as 



Q+ {X)g+ (k) ^^^^^^ y^ UkUk' hk (X)hk' (k) 



C(A,k) 



M 



{lifiKf 



k,k' 



4Me' 



^g4A)g.M+ '^"''\'+'^:-w .(58) 



A + K 



with the equilibrium correlation function C^'"^(A) in the 
Laplace domain defined in Eq. (|C12p . In App. [C] it is 
shown that hk{X) can be written in terms of Q^ (A) and 
^+^(|i^fe|) [see Eq. (jClip j. For uj = ujo wc use the fact 
that Uk = g+{\vk\) and we find from the definition of the 
equilibrium correlator 



''-^'>-wV^ 



A2 



g+(A)-^+(|z.fe|) ,(60) 



which is derived in App. [C] [see Eq. ()C12|) ]. that the 
desired non equilibrium correlation function of a quan- 
tum Brownian particle with initial position measurement 
reads 



A/f A' — A ~ 
C{X,n)^-—^C''^'^{X)C'^^{n) + 



h A2 
ci''q(A)-hCi°q(K) 

X + K 



^^5^iX)g^{n) 



(61) 



The classical correlator for an initial position measure- 
ment with outcome g*^ = and uncertainty A can be 
obtained from Eq. (|B8p by replacing (q^ ) by A^. The 
initial momentum is not measured and it is therefore dis- 
tributed according to the Boltzmann-law with (uq) = 
(/3M)~^. In the limit of a sharp position measurement 
(A — > 0) the classical correlator of an equilibrium particle 
reads C^^X, k) - ^Muj^ Ci°q(A)Ci'=i(K). As to the quan- 
tum correlator we note that A' = for A -> and the 
sum of the first and the third term in the rhs of Eq. (pT|) 
yields C''q(A,K) - M/{hA) C'^'"i{X)C^'"i{K) which already 
has the form of its classical counterpart. It is easy to show 
that in the high temperature limit /3fi <C \uj'^ — 7^/4| the 
two expressions coincide exactly. The role of the second 
term in the rhs of Eq. (|6T|) remains to be discussed: it 
describes the initial momentum due to Heisenberg's un- 
certainty relation. Consequently, it diverges when the 
initial position measurement becomes sharp unless one 
considers that fi/^Me^) = h^/iM^A"^) -^ even though 
A -> 0. More precisely, when f3h\uj'^ - 7^/41 ^/^ < 1 
we have C^^^ ^ 1/{(3Muj'^) and 5+ ^ 1/uj. In order for 
the second term of the rhs of Eq. (|CT|) to be small com- 
pared to the other terms the condition A ^ At must 



hold, with At = J [ifi^ /{2-kM) the thermal de Broglie- 
wavelength of the particle. Only then can one speak of 
a classical particle: the condition j3fi\Lo'^ — 7^/41 -^/^ ^ 1 
that properly defines the "high temperature regime" is 
not sufficient. One also has to take the "macroscopic 
measurement limit" defined through A 3> At- 
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In the real time domain the quantum correlation func- 
tion is 

C{t,t') = ^^^C^^^{t')C'''^{t') 

+ ^g+{t)Q+{t')+C^^%\t-t'\), (62) 

where C^°'^{t) is the real time equilibrium correlation 
function. This result applies to any kind of spectral den- 
sity of the bath. We will use the correlation function 
Eq. (|62p in Sec.|V]to study the non equilibrium dynam- 
ics of an impurity in a Luttinger liquid bath for which a 
specific spectral density of the bath applies. 

We are interested in the equilibration behavior of 
Eq. (|62p in the presence of Ohmic dissipation for which 
the spectral density behaves as S{ui) ~ 70; for small 
w, with 7 playing the role of a friction coefficient [for 
more details see App.ID]. It is of special interest to study 
the strong quantum regime (ifi ^ |w^ — 7^/4|~^/^. By 
using the long-time limits Eq. (|D5|) and Eq. (|D6|) pre- 
sented in App. |D] valid for t ^ 7^^ we find that the 
propagator Q^ (t) exponentially approaches zero whereas 
the equilibrium correlation function C^'^'^{t) relaxes with 
a power law r^ t~'^ for i — >■ 00. In the long-time 
limit the correlator Eq. (p^ thus relaxes as fast as the 
squared equilibrium correlation function. Therefore, for 

/3?i>|Lj2_^2/4|-l/2 . 



domain reads 



cit,t') = c^'^^ilt - 1'\) + o [{tt'y 



(63) 



when t,t' ^ 7^^, and the equilibrium function C^°'^{\t — 
t'\) = C°^{t,t') is asymptotically approached during the 
algebraic relaxation of the non-equilibrium terms. Con- 
sequently, care has to be taken in experiments when an 
equilibrium system is desired at very low temperatures af- 
ter an initial position measurement. The relaxation of the 
system is slow independently of the dissipation strength 
7. At high temperatures pti <^ jw^ — 7^/4|^^/^ the relax- 
ation is of order C'(e~'^*) [see the discussion in App. |D] 
and therefore exponential as in the classical limit. 



B. Quantum quench in the confining potential. 

In this section we desire to gain insight into the equi- 
libration process of a quantum Brownian particle after 
an abrupt change in the trapping frequency. At i < 
the particle is confined in a harmonic potential with fre- 
quency wq > 0. At t = the experimentalist abruptly 
changes the strength of the harmonic potential resulting 
in a higher or lower trapping frequency. We do not con- 
sider an initial position measurement since we assume 
that the particle is already localized by the initial har- 
monic trap. Hence we set A = A', H. = fl' , Xm = Xm = 
and e ^ 00. By starting from Eq. ((55)) and by using 
Eq. PO]) it is straightforward to show with the methods 
employed in App. [C] that the correlator in the Laplace 



C(A, k) 



A/gO(A)^0(K)A + K 






(A) 



^loq 



in) 



(64) 



where 5+(A) = l/[A^-|-A7(A)+Wo] is the propagator with 
initial frequency wq, and 



Oa) = ^E' 



1 

(3AI 



A2 



gl{x)-gUWk\) (65) 



is the equilibrium correlation function of a particle in a 
harmonic potential with frequency luq [see Eq. (|C12[) ]. 
The structure of Eq. (|64l) is very different from the clas- 
sical counterpart Eq. (|B7|) . Still, by using the high tem- 
perature approximation Cq°'^{X) ~ — ^+(A) — 1/ujq /A 

[see the end of the discussion in App.[D] one recovers the 
classical expression Eq. (jB3p . 

The equilibration time for Ohmic dissipation in the 
strong quantum regime /3h 3> \uj^ — ^^/'i\~^^'^ can be 
found in the following way. Note first that when Cq°'^{X) 
is multiplied by G+iX)/Q'^{X) one obtains a new function 
that we call C: 

X 



'''''-w^< 



A2 



G+{x)g+{\vk\) 



(66) 



Second, we observe that the Laplace transform of [dt + 
dt')f{t,t') is equal to (A + k)/(A, k) - f{t = 0; k) - 
/(A; t' = 0) where /(A, k) is the Laplace transform of 
the generic function f{t,t') with respect to both t and t' 
while /(A; t' = 0) [f{t = 0; k)] is the Laplace transform 
of /(i,t' = 0) [f{t = 0,t')] with respect to t (f). Now, 
by choosing f{t,t') = C{t,t') we find for the correlation 
function Eq. (|64p in the time domain 

(dt + dt')cit,t') = c'it) [1 + iuj' - iul)g+{t')] 
+c'it')[i + {oj'-Lol)g+{t)] 

-C(t,0)-C(0,i') • (67) 



From Eq. (|64p we easily find the expression of the Laplace 
transform of C(t, t' ~ 0) by multiplying C(A, k) by k and 
by taking the limit /c — > 00: 



C(A; t' = 0)= Coi°'i(A) [1 + (c.2 - ujl)g+iX) 
= C'(A). 
Accordingly, Eq. (|57|) simplifies to 

{dt + dt>)cit,t') = c'(t)(c.2 - i,l)g+{t') 

+C'(t')(c.2-c.2)g+(t) 



(68) 



(69) 



Finally, it is clear that {dt+dt')C{t, t') = when C(t, t') = 
C'''i{t,t') = Ci'='i([t - t'\) so that the terms of the rhs of 
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Eq. ((69|) can be understood as the derivative of the non 
equilibrium part of C{t,t'). In combination with the re- 
sults of the discussion in App.|D]the equilibration behav- 
ior of the correlator at low temperatures after a quench in 
the trapping potential can be summarized as follows [by 
noting that the asymptotic long-time behavior of C'{t) 
and C^'"^(t) are identical]: 



C(t,i')=Ci°'i(|t-t'|) + 



-'It/2 



A' 



-7*72 



/t 



(70) 



for t,t' ^ 7~^- Consequently, the non equilibrium con- 
tributions are exponentially suppressed which leads to a 
faster equilibration than the one found after an initial 
position measurement. 

At high temperatures the relaxation of C'{t) and 
C^^^^it) are both exponential of 0(e~'^*/2) gQ ^^^^ -^^ 
the high temperature regime C{t,t') = C^'"i{\t - t'\) + 



V. EVOLUTION OF AN IMPURITY IN A ID 
ROSE GAS 

Many ID systems are described by the Tomonaga- 
Luttinger theory J37l - l39l |. The evolution of impurities 
in such a Luttinger liquid (LL) has attracted much at- 
tention in recent years |17l . |l9[ since this problem can be 
explicitly realized with cold atom systems. In particular, 
modern techniques allow one to tune the interspecies in- 
teraction strength J4(>|42| | so that it has become possible 
to study the diffusion of a minority species within an en- 
semble of majority atoms, as a function of the interaction 
and the trapping potential [l6| . 

In this part of our work we apply the non equilibrium 
formalism developed in the first part of our article to 
such an impurity-LL system. In particular, we seek to 
mimic the experimental process described in [ig with our 
theoretical description. In this experiment the impurity 
atom is trapped in a ID harmonic potential together with 
an ensemble of a different kind of atoms that form the 
LL. The impurity is initially localized at the center of 
the confining potential by a laser blade. When the whole 
impurity-LL system reaches equilibrium the impurity is 
released. The equal-times position correlation function 
of the impurity, C{t,t), then shows damped oscillations 
which strongly suggest that the impurity is de facto a 
quantum Brownian particle moving in a quantum liquid 
bath. 

In the following we will present a precise description of 
the impurity motion in the LL from the quantum Brown- 
ian motion point of view. The LL itself will play the role 
of an exotic quantum bath that we here characterize. 



A. The impurity model 

The impurity and the atoms constituting the bath are 
all confined in a harmonic potential. We therefore take 



the Hamiltonian of the impurity, Hs, to be of the stan- 
dard form Eq. ([1]) without external force {H = 0) and 

with the harmonic potential V{q) = ^^^cf'. We assume 
that the interaction Hamiltonian between the position 
operator q of the impurity and the density of the boson 
liquid is of the form 



n 



SB 



dxdy U{x - y)p{y)6{x - q) , (71) 



with the density operator p{x) of the LL approximately 
described by 



p{x) ~ po{x) V0(x) , 



(72) 



where po{x) is the unperturbed density of the fiuid in the 
ID trap and 0(x) is the density variation. The Hamilto- 
nian of the free LL reads 



n. 



h 

2^ 



dx 



uKr 



{T,tl{x)f 



Kl 



{vmf 



(73) 



where H(x) and 4'{x) are conjugate operator fields. Equa- 
tion ([75]) describes the low-energy properties of a Lieb- 
Liniger gas |43l | with a contact interaction potential 
hvLS{x). The parameters u (with the dimension of a 
velocity) and Kl (dimensionless) have to be determined 
numerically for general vl- In order to reduce the com- 
plexity of the problem we will assume the background 
density po to be constant in the following. Accordingly, 
we define 



1 



Po 



dx po{x) , 



(74) 



where L is a length scale of the order of the length of the 
trap. Note that in this modeling we have not added the 
quadratic confining potential to the LL. 

Since the Bose gas is confined in a space of length L, 
the wave vectors are quantized with values fc„ = im/ L 
with n an integer. The Fourier representation of Eq. (j7ip 



n 



SB 



^E^- 



pik„q 



Vl 



" In 

(P[Kn 



(75) 



where we used S{x — x') = (1/i) X^n ^''^"^^"'^ "* ^'^'^ ^^ 
neglected a constant contribution. We assume that the 
potential U has the form 



Uk = five' 



fc|/fce 



(76) 



with some finite cutoff kc that depends on the micro- 
scopic properties of the interaction. The parameter v 
has the dimension of a velocity and it determines the 
strength of the impurity-bath interaction potential. 
After redefining the fields according to 



and 



[x) ^ V{nKL/n) i{x) 



ti{x) ^ Vih/nKL) n(a 



(77) 



(78) 
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one introduces the bosonic ladder operators 



(79) 



which describe bosonic density wave excitations with 
sound velocity u 4J]. The full Hamiltonian now takes 



the form of the FrohHch polaron Hamiltonian, which in 
the second quantization language reads 



and 



^--Vi!(^(^)-^^(.)), 



(80) 



fce{fe„} 



Moj 






2Tl\k\ 



1/2 



Uk 



[ike''"'') bl + [ike''"")* b. 



(81) 



r 



For each k mode the coupling between the operator e**''^ 
and the bath operators b\. and 6fc is bilinear, so we can 
use the general results derived in Sec. IIIIBI bv consider- 



ing that e**''^ represents a different operator for each k 
that is coupled to one harmonic oscillator. The resulting 
influence functional is hence a sum of one-particle-one- 
oscillator influence functionals. By using Eq. (|A27p we 
thus find 



Hq^^q-.q""] 



Y. { ^'^ Aa Kk{-iT + ia)e'^''"^''^e-'^'i^'''^ 



dr / ds Klkis-iT)e' 





— / ds I du 

/o Jq 



^~ikq+(s) _ g-ifcg [s] 



^ikq+{s} _ gifcg (s) 



K^k{s - u)e 



—ikq^ (u) 



~K\{s-u)e 



-ikq [u] 



■ (82) 



r 



Note that, since in the initial Hamiltonian there is no 
counter-term the counter-terms in Eq. ()A27|) are absent 
in Eq. ([5^ . The fc-dependent kernel reads 

1 /u|fc| ~ 2 cosh[u\k\j3h/2 ~ ieu\k\] 

I JJu I 



Ku{e) 



(83) 



27rn ' ' sinh[M|fc|;3?i/2] 

A well known feature of the Frohlich polaron Hamil- 



tonian (|5l]) is its polaronic mass shift which leads to an 
effective impurity mass M* > M greater than the "bare 
mass" [l[. Another process described by this Hamilto- 
nian is that during a collision between the impurity and 
an LL atom the former loses momentum fik by creating 
a density wave excitation b\, in the LL. However, the LL 
is itself confined in a harmonic potential and one may 
expect to have some momentum transfer absorbed (or 
provided) from the LL to the optical trap. We will here 
simply assume that the spring constant of the optical 
potential is renormalized by such a process in such a 
way that it balances the polaronic mass shift and we will 
henceforth work with the bare Hamiltonian ((5T|) . 

If the oscillations are small (if the impurity potential 
is sufficiently steep) we can expand the e^^'^ in Eq. 



to second order in k. Note that the linear order vanishes 
in Eq. (j82[) due to the symmetry k t-^ ~k. Ignoring the 
zero-th order in k we can make the replacements 



g»fcq«(r)g-jfcgO(^) ^ 



[q"{r)-q'{a)]\ 



^lkq\T) 



^~ikq+{s) _ g-ikq (s) 



H> 



fcV(r) [qH^)^q-{-^)\--:^[q^{sf-q-{sf] , 



(84) 



e [q+{s) - q-{s)\ [Kk{s - u)q+{u) - Kl{s - u)q- {u)] 
-^[q-'{sf-q-{sf]{Ku{s~u)-Kl{s~u)) . 



13 



Steep potential: The Luttinger bath and the 
initial condition. 



If the external potential is steep we expect the Gaus- 
sian approximation of the Frohlich Hamiltonian ([5^ to 
be valid. Up to this order in k eaeh fc-mode plays the role 
of one bath harmonic oscillator. Accordingly, for i — > oo 
we define the spectral density as 



^H = 7 E 



KM' 



fee{fc„} 



2TTfl 



\lJu\^ 8{u, - u\k\) 



TT/i / LO 



(85) 



with LOc = u\kc\/2 and 



AHKlv'uj 



2, ,3 



(86) 



Hence, the bath induces super-Ohmic dissipation with a 



power law behavior S{lij) 



for small 



This is the 



main difference from the analysis presented in Sec. IIIIl 
In terms of the fundamental kernel [e.f. Eq. (|A28p ] 



Ki0)= I d^5(^)POshK/?V2-^e)] 



sinh[ujl3h/2] 



(87) 



ph 



and with /^ dTK*{s 



= AI'-f{s) and 



/p du [K{s - m) - K*{s - u)] = -iM-f{0) + iMj{s), 
where 7(5) is defined in Eq. (|A42[) . we obtain 



(•/3R 



1 rPiT- „ rpn- r^ 

$[9+,(?-,q°] = - / drda K{-iT + ia)[q"{T)-q°{a)] -I dr ds K* {s - tT)q° (t) [q+ (s) - q~ (s)] 

4 Jo Jo Jq 

^^■^7(0) r, r +. ^2 -, ^2l 



ds / du [q+{s) - q-{s)] [K{s - u)q+{u) - K*{s ~ u)q- {u)] + 
Jo 



dslq+isy-q-isY 



The first line in Eq. (|55)) shows that only the non-local 
part of K{—iT + ia) contributes. Hence we can re- 
place the first line by ^ /^ drdcr k{—iT + ia)q°{T)q°{a). 
Furthermore, we see that the last line in Eq. ([SS]) ex- 



I 

actly represents the counter-term proportional to /.j. By 
rewriting Eq. (|55)) in terms of the kernels Kfj and 7 [see 
Eqs. (|A33[) and (|A42[) ] and using the transformed vari- 
ables 2; = {q~^+q~)/2 and x = q^ —q^ the action becomes 



J 



./3 a 

^[x,x,q°,Xt,Xf,Xt]=i dr 

Jo 



M.Q, .2 MUJ^ „. .2 1 '^^''■ 

y9°M +^-g"(r) +- 



^9"(r)' + '-^q"{rr + ^ I da fc(r - a)<z°(r)g"(a) 



Jo 

t ft 



+ : 



ds / du Kji{s — u)x{s)x{u) + 



ri3h r-t .t 

+ I dr I ds K*{s - iT)q"{T)x{s) + / ds 

ih 
2A2 



AIx{s)x{s) — Mu x{s)x{s) ~ Mx{s) — / du 7(s — u)x{u) 

ds Jo 



(89) 



We remind the reader of the main approximations used 
so far. First, we used the Gaussian approximation ([84|) 
of the Frohlich Hamiltonian. Second, in order for the 
action ([55]) to make sense we assumed ut to be large 
enough, so that non Gaussian effects are not too impor- 
tant. Note, that the mass and the potential renormaliza- 
tion can modify the oscillation amplitude and the final 
width of the impurity position. Finally, we interpreted 
the laser blade that initially localizes the impurity at the 
center of the quantum liquid as an initial position mea- 
surement with outcome <?,„ = and uncertainty A. The 
effect of the initial position measurement is incorporated 
into the action via the last term of the rhs of Eq. (|89|) . If 
the localization is performed itself by a very steep trap- 
ping potential with frequency flo the particle is in its 



I 

harmonic oscillator (with respect to flo) ground state at 
initial time. We then have 



= 2f7o 



(90) 



This approximation is disputable since one could also 
consider the initial localization of the impurity as stem- 
ming from an initial trapping potential. The subsequent 
release of the impurity would then rather be described by 
a quench in the harmonic potential [see Sec lIVB] . How- 
ever, in real experiments the "high temperature" regime 
fi/3u! ^ 1 prevails. From the discussion in Sees. IIV Al 
and llVBl we know that in this case the difference between 
the particle motion after an initial position measurement 
and the one after an initial localization due to an initial 
trapping potential (followed by a quench in the potential) 
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is blurred. Since it is technically easier to deal with a po- 
sition measurement we prefer this method to a quench 
in the trapping potential. Note that the "high temper- 
ature" regime is not equivalent to the classical regime, 
as we pointed out in Sec. IIV Al as it does not fulfill the 
"macroscopic measurement" condition A ^ At with At 
the thermal de Broglie- wavelength of the impurity. For 
more details go back to the discussion in Sec. IIVAI 



C. Signature of a Luttinger liquid bath 

A typical experimental scenario consists in holding the 
impurity in the center of the trap at i < and switching 
off the localizing potential at t = to let the impurity 
move in the residual harmonic potential. Let us for a 
moment forget about the potential renormalization and 
the mass shift. By using the results found in Sec. IIV Al we 
find in the limit e — >■ (which corresponds to an almost 
perfect localization of the impurity at t = 0) 



C(A, k) 



Ule 



M ~ 

^g+{\)g+{K)-—c'^^{x)c''^'^{n) 



+ 



1 



A + K 



hA 



(91) 



In the time domain the correlation function reads 
Cit,t') 



^G+m+it') - ^c'^'^itr'^'^it') 



-Ci°i(|i-f|) 



(92) 



For the moment, experimental measurements focus on 
the equal-time correlation (which corresponds to the 
time-dependent variance of the position) 



C(f,t) 



4A/e2 



5Ut) 



M 
hA 



C^^^^itf + C^'='i{0) . (93) 



The formula Eq. (|93p is valid for all kinds of baths and 
for very small polaronic effects and potential renormal- 
ization. 

In order to gain further insight into the dynamics of 
the impurity we need to understand the contribution of 
each of the three terms in Eq. (|93|. 

We start from the analysis of the propagator 0+ (t) . To 
determine its time-dependence we need the specific form 
of the spectral density of the LL bath (|55|) . In terms of 
the function 



K^) 



1 

v2 ^4 



dC C' e 



< 



z 

y 



z- 

y 



dC 



C 



e + z 



C^+Z2 



(94) 



the Laplace transform of the damping kernel 7(A) [see 
Eqs. (jA42p ] can be recast as 



A7(A) - ^5(AK) 



(95) 



and the propagator (which is proportional to the linear 
response function) reads 



Q+W 



1 



A2 + i^g(AK)+c.2 



(96) 



Our objective is to find the oscillation frequency and the 
damping of the impurity motion in the small coupling 
limit. Thus, we need the inverse Laplace transform of 
Eq. (j96p which can be expressed in terms of the Bromwich 
integral. 



^+(0 



1 

27ri 



C+'icxD 



dA e^* ^+(A) 



(97) 



where the real number c is greater than the real part of all 
poles of Qj^(X). The integral in Eq. (|97p can be solved by 
displacing the complex contour towards the left and by 
evaluating the encountered residues. Hence, we seek for 
the complex points at which the denominator of Eq. ([96|) 
vanishes, that is z^ -I- fi'g{z) + w^/a;2 = with 

z = i(T — T, where — <C 1 and n' = 5- . (98) 

Note that z, cr, F and fi' are dimensionlcss and a and F 
real. From Eq. (^7)) it is clear that F corresponds to the 
exponential damping (measured in units of Wc) of the 
propagator and hence F has to be positive (whereas a 
is the oscillating frequency measured in units of Wc and 
can be positive or negative). Thus, the first singularities 
encountered when displacing the contour towards the left 
are the imaginary axis for which g{z) is obviously singular 
when Re z = 0. In other words, we have to account for 
the difference g{ia + 0+) — g{ia — 0+) when displacing 
the integration contour past the imaginary axis. We will 
come back to this point in just a few lines. The integral 
in the second term of the rhs of Eq. (IMl) can be recast as 



/(-F,a)^ 



(,^ -I- z^ 



dC 



-kIC 





I'^l ./o "^^ C' - 1 - 2iF/a 
"F" 



1 



A{\a\)+ sign 



IFI 



-(l + |a|)e- 



a 



kl 



!!L,-ki 



I sign 



s(H) 



where we expanded the denominator of the integrand to 
first order in F/cr and we defined 

A{a) = sinh a Chi a — cosh a Shi a , 
B[a) = a [cosh a Chi a — sinh a Shi a\ 
+ cosh a Shi a — sinh a Chi a , 

where Shi and Chi are the hyperbolic sine and co- 
sine integrals, respectively. The contribution from the 
crossing of the imaginary axis is now easily obtained: 
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A/ = l(0+,cr) - I(-0+,cr) = -11X6-^"^/ a. After adding 
—z'^/S.I/2 to g{z) we thus obtain for Rez < 



9{z) ^ - 



^ iTTO^ 1^1 a 



2 

7rra2 






(99) 



where we neglected all real terms of 0((r/cr)^) and all 
imaginary terms of 0{T/(j). 

At the poles, the real and the imaginary part of the 
denominator z^ + ^' g[z) + {ui / uicY vanish simultaneously. 
By using the form of g{z) just derived we see that this 
can be achieved with the choices 






(100) 



- ^ ^2 ^ /^ r^3^(^) + ^2] ^ 2r2(5 + |a|) . (101) 



Note the symmetry a i— > —a in Eq. (jlOip . Equa- 
tions (jlOOp and (|10ip determine the oscillating behavior 
of Q+{t): r corresponds to the damping and a to the 
frequency of the oscillation (both quantities measured in 
units of Wc). We assumed the damping F to be small 
compared to a which by virtue of Eq. (|100p translates 
into 



4 4e 

H' < — e" > — w 3.461 

■KG TT 



(102) 



The behavior of the oscillating frequency a (measured 
in units of Wc) strongly depends on the trapping fre- 
quency uj/lUc- Figure [T] shows the dependence of |ct| on 
fi' = fi/{Muj'^) for the cases uj/uJc = 1.0 and uj/uJc = 3.0. 
When h(3uj <C 1 (as in [1^]) the properties of the oscil- 
lations of G+it) coincide with the ones for C^°'^{t) [see 
App. [B]. Then, by observing Eq. (|93)) one deduces that 
the actual oscihation frequency oi C{t,t) and its damp- 
ing are rather 2a and 2F, respectively. As can be seen 
for moderate to high trapping frequencies there is an in- 
crease of the oscillator frequency followed by a decrease 
induced by the bath [see Fig. [T] . Such a behavior is not 
observed for small frequencies (compared to ujc). The 
experimental results in [la ] confirm such a "peak" al- 
though very few data-points are shown in this paper and 
the error bars are quite large so it is hard to draw firm 
conclusions on the actual experimental behavior at this 
stage. The "peak" in the curve of the oscillation fre- 
quency can be used in further experiments to determine 
whether the bath is actually described by a LL or not 
since it is a direct consequence of the non-Ohmic spectral 
density (|85l) . In Fig. [2] we show the correlator Eq. ((93)) 
for values of the impurity-bath coupling (/x' = 0.5 and 
fi' — 0.2) in the "high temperature" regime /3?iw <C 1 
typical for experiments. The parameters were chosen to 
be l/(4e^ajc) = 10 in both images for the thick lines, 
Lo/ujc = 1.0 in the upper and ui/coc = 3.0 in the lower im- 
age. The thin line in the upper image has been obtained 
with l/(4e^ajc) = 50. The curves qualitatively agree with 



3.02 



3.00 



2.98 




FIG. 1. Dependence of the oscillator frequency a on the 
impurity-bath coupling fi' — jj./{Mlj'^) for lj/ujc = 3.0 ob- 
tained from Eq. (|101|l . A logarithmic scale has been used. 



the experimental data [Iq . One clearly recognizes the in- 
crease of the oscillation frequency induced by an increase 
in the coupling fi' for uj/uJc = 3.0. The typical oscillation 
width in experiments is about 15//m with a frequency of 
a ~ 550s~^. The final position width of a (high temper- 
ature) quantum Brownian particle is given by l/(/?Mcj^) 
[3]. With the mass M of the '^^K atoms used as an im- 
purity and the bath temperature T w 350nK in [l6[ this 
yields ^jC{t,t) K, 15/xm which is even quantitatively the 
right value. 



VI. SUMMARY AND OUTLOOK 

Motivated by recent experiments on impurity dynam- 
ics in quantum liquids we studied the non equilibrium 
dynamics of a quantum Brownian particle coupled to a 
quantum thermal bath of harmonic oscillators for generic 
Gaussian initial conditions. We derived a closed expres- 
sion for the non equilibrium correlation function, which 
we showed to be easy to derive as variations of a gen- 
erating functional Eq. ([5^ . We used the analysis in Q 
as a starting point to obtain this generating functional 
by employing path integral methods. We then showed 
that factorizing initial conditions (where the bath and 
the particle are initially uncoupled) are a special case of 
the non-factorizing initial conditions on which a position 
measurement has been performed. We demonstrated the 
correctness of our approach by deriving the equilibrium 
correlation function without imposing time-translational 
invariance (presented in App. ICj) . 

We applied this general formalism to the study of three 
physical situations. First, we studied the equilibration 
process of a trapped particle after an initial position mea- 
surement. In this case we considered Ohmic dissipation. 
While the classical (high temperature) correlator relaxes 
exponentially on a time scale 7~^, the low temperature 
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Moj^C(.t,t)/h 
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Miji)cC(t,t)lh 
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FIG. 2. Theoretical results for the variance of the impurity 
position C{t,t) after its quantum Brownian motion induced 
by a LL bath in the "high temperature" -regime which is 
experimentally available. We chose l/(4e^a;c) — 10 for the 
thick lines in both images, and lo/ujc = 1.0 in the upper and 
Lo/uJc = 3.0 in the lower image. The impurity-quantum liq- 
uid couplings are ^' = 0.5 (thick straight lines) and ^! — 0.2 
(thick dotted lines), respectively. In the upper image the thin 
line has been obtained with l/{4:e'^iJc) = 50 and fi' — 0.5. 



correlator C{t,t') in the strongly quantum regime shows 
an algebraic relaxation of the form l/(ti')^ which is in- 
dependent of the dissipation strength 7. Therefore, the 
information that an initial measurement on the system 
has been performed persists for a very long time. We 
then showed that the equilibration process is different if, 
instead of a position measurement, a sudden quench in 
the trapping potential is performed at the initial time. 
We showed that in this case the relaxation is exponential 
in time even in the quantum regime with the slight differ- 
ence that at very low temperatures the relaxation time 
is of order 27"""^ rather than 7^^ for high temperatures. 
Accordingly, the relaxation due to quantum fluctuations 
is almost as effective as thermal relaxation in this case. 

Third, we applied our approach to the motion of im- 
purities in ID ultra cold quantum gases confined in a 
very elongated harmonic trap. For low excitation en- 
ergies the ID gas is well described by the Luttingcr- 
theory which predicts a Gaussian theory. We modeled 
the impurity-gas interaction by a simple contact poten- 
tial. The resulting Hamiltonian is also known as the 
Frohlich polaron Hamiltonian which features a non-linear 
coupling of the polaron (here the impurity) to the bath 



(here the Luttinger liquid). It is well known that this 
Hamiltonian induces a mass shift on the impurity de- 
pending on the interaction strength. Although the mass 
shift is clearly observed, no significant change of the os- 
cillation frequency is detected in [ig. Hence it is likely 
that the confining potential is renormalizcd as well in 
such a way that it counteracts the polaronic mass shift. 
We then proposed to use, as a first approximation to 
this problem, the Gaussian approximation of the polaron 
Hamiltonian (more subtle issues linked to polaron and 
harmonic trap renormalization together with a quantita- 
tive comparison between theory and experiment will be 
presented elsewhere) . The spectral density of the result- 
ing exotic "bath" turns out to be strongly super-Ohmic 
with the low frequency spectral density S{u!) ^ uj^. With 
the help of the generating functional approach we were 
able to deduce the non equilibrium correlation function of 
the impurity position. Our formula correctly reproduces 
the oscillating behavior of the position variance (which is 
the equal time correlation function) as well as the initial 
momentum of the impurity due to its strong localiza- 
tion at the center of the quantum liquid which is a pure 
quantum effect. We showed that the typical impurity 
position width measured in [16| quantitatively matches 
the theoretical value for a quantum Brownian particle. 
This suggests that non equilibrium quantum Brownian 
motion theory can describe many aspects of such impu- 
rity dynamics. The super-Ohmic dissipation leads to a 
curious phenomenon: upon increasing the coupling be- 
tween the impurity and "bath" the oscillation frequency 
increases before decreasing. This peculiar behavior is in 
strong contrast with that induced by Ohmic dissipation, 
for which the damped oscillator has a lower frequency 
than the free one and it can be used as a signature for 
Luttinger liquid baths. 

To our knowledge our formalism is the only one that 
can fully access non equilibrium correlation functions in 
quantum Brownian motion in contrast to the density ma- 
trix approach in [9[ that yields only one-time quantities. 
It proved effective for the description of impurity dynam- 
ics in Luttinger liquids. By using a quantum Langevin 
equation the authors of [1^ have described the impurity 
motion with some success. However, the same authors 
pointed out that the Langevin approach failed to grasp 
all aspects of the experimental evidence. We are con- 
vinced that one can fit our theoretical results to experi- 
mental data as soon as more data become available in the 
future. The parameter ujc which depends on the actual 
impurity-quantum liquid coupling can then be experi- 
mentally determined within our framework. We hope to 
address this topic in a future study. 

To conclude, our theoretical framework proved effec- 
tive to deduce all kinds of non equilibrium correlation 
functions and we obtained promising results for its ap- 
plication to impurity motion in ID systems. Further pos- 
sible applications of our approach include various topics 
such as the study of decohcrence problems, the impu- 
rity motion in other types of systems or non equilibrium 
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quantum transport phenomena. The present work is ex- 
pected to provide a further motivation for a thorough 
understanding of quantum Brownian motion. 



Appendix A: Coherent state path integral 
formulation 



In this Appendix we detail the derivation of a path- 
integral representation of the generating functional of the 
multi-time correlation functions of a quantum particle in 
contact with a generic quantum bath made of an ensem- 
ble of harmonic oscillators. 

The terms in the Hamiltonian of the dissipative quan- 
tum Brownian particle that depend on the bath variables, 
Eqs. ([2]) and ([3]), can be rewritten in terms of creation 
and annihilation operators of the bath oscillators in view 
of a later use of coherent states. One defines the creator 
and the annihilator of the n-th oscillator mode by 



^t 



The operators d„ and a+ satisfy the bosonic commuta- 



IrrinUJri 


i 


(Al) 


2h 


Xn Pn 
mnOJn 

i 


IrrinUJn 


(A2) 


2h 


Xn 1 Pn 
mnOJn 



tion relations [a^ , i 



and [a+ , d. 



+ ;;+! 



[diit ^rn\ 



0. Equations ^ and ^ read in terms of the d^ and dn 



Hb 


= ') ' ?iw„d+d„ , 

n=l 


(A3) 


-isB 


= ^ 5n9(an + On) ■ 


(A4) 



Here we introduced the notation g„ = ^/hc^j2mnOJn- 

We introduce the coherent states of the harmonic oscil- 
lators, which are particularly suitable when dealing with 
the bosonic ladder operators in Eq. (|A1[) . 
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|0), 



(CI = (0| 



.Cs. 



(A5) 



where ^ is a complex number and ^* its complex conju- 
gate, d"*" and d stand for the creation and annihilation 
operator of each harmonic oscillator. The coherent states 
are eigenstates of the annihilation operator, that is 



(cir 



(A6) 



with the properties 

(e|C)=e«''; and 1' = J dCd^e-^^^lOi^A^) 

where 1' denotes the unity matrix of one oscillator. 
Hence, the unity matrix of the whole particle-bath sys- 
tem 1 can be written in terms of the product states 
k, {?«}) as 



The trace of any observable 13 that depends on the par- 
ticle and the bath operators can be expressed as 

TtB = [dq fH {dCtdCn e-«"«" } {q, {^n}\B\q, {6.}) ■ 

(A9) 
The generating functional can now be obtained by sup- 
plementing the potential V{q] s) in Eq. ([T]) by a linear 
term —H{s)q where H{s) is a c-number function that 
plays the role of an external source. To be more explicit, 
we introduce two distinct sources H{s) and H'[s) in the 
potential V for the time evolution operator 

l<^{qf,{Cn,f};qi,{inA''t) = 

('?/,{W}|te-'^/o'^*'«(*')|g„{e„,0) 

(AlO) 

and its Hermitian conjugate 

(9UC,}|ttenj^*'«(0|g^,{^;_^}), 

(All) 

respectively. T"^ is the anti-chronological time order- 
ing operator. In the following we will use the short- 
hand notation dfi = O^i s~^"''^"''d^n A^n,i ■ All cor- 
relation functions of the position q can be obtained by 
taking the corresponding variations of the trace of the 
time- dependent density matrix p{t) = A^po/C*, Eq. (|16p . 
with respect to H{s) and H'{s). The matrix elements of 
p{t) = ICpolC* are given by 

w{qf,{^nj};q'fAC.f};t) - {qfA^nj}\ m k/,{Cj}) , 

(A12) 

^^iq^,{U^hql{CA) = (q^dCnAl po HACA) ■ 

(A13) 

The path integral representation of /C and /C* are found 
by using standard methods: 



^qf,{^nJhq^,{in,^}■.t) = J^^q^^^i 



xexp(-5[g+,{e}] 



ic* iq'fACnjh qi {CA; t') = J 'Dq-vc 



(A14) 



xexp( --5*[g ,{C„}; 

(A15) 

where we make clear by the superscripts + and ~ which 
path belongs to K. and which to /C*, respectively. The 
real time interval [0, t] has been discretized into T G IN^ 
steps of length At with t = AiT. The functional in- 
(A8) tegration measures are defined as Vq = HiJi dqj with 
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qj = q{jt/T). The terms contributing to the total ac- 
tion, S[q,^ n] = Ss [q] + Sb [Cn] + SsB [g, Cn] introduced in 
Eqs. ()A14[) and (jAlSp . read in discretized form 



Ss[q]=Y.At 



M fqj-qj^i 



-HijAt)q,] , 



At 



T-l 



^ -V{q,;jAt) 



(A16) 



1. Integration over the bath variables 

The influence functional in discretized form reads 



^si^,] +5sB[g,C„] = ih^£,lj{^nj -£.n,j-l) 



n j = l 

X e 



igr. El^T+i At'g,(4„,,+C:,,)-^g„ EflsT+i Atg, (?„,,+«;,,) 



i=i 



(A21) 



(A17) 



The reduced density matrix depends only on the par- 
ticle variables and the external sources, 



W{qf;qy,t) = J dq,dq'Mde^d^f JC{qf A^njh Qr, {^nA^t) 

X W(gi,{C„,i};g,-,U^,J)/C*(gy,{C„,/};g-,{^'j_J;t) . 

(A18) 



which depends on the correlation matrix 
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/I ••• 

-fc2 1 ••• 

-fca 1 ••• 

-ki 1 
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-ksT 1 / 



(A22) 



In Eqs. (jA14p and (JAISP we omitted normalization fac- 
tors that do not depend on the bath nor on the particle 
variables. Note that the integral runs only over bath and 
particle variables with an index between 1 and T—1 since 
Qo — li; Qt ~ <lf (and analogously for the bath variables) 
arc fixed for K. and qo = q'f, qr = q'^ (and analogously for 
the bath variables) are fixed for /C*. 

The path integral description of po is obtained by di- 
viding the imaginary time interval [0,/3?i] into T time 
steps. Consequently, by using p = e~°'^° we find 



W{qi{en,^hq^AU^}) = 

Jvq'^Yl^VCe.p^-^-Soiq'^C] 



1 + iAtuj for j <T, kj =1- iAtoj for j > 2T 
and fcj = 1 - At'w for T < j < 2T. N ote th at in this 
notation there are ST time steps in Eq. (|A21|) and qj = 
q+ for j < T, qj = g^" for T < j < 2T and q^ = qJ for 
j > 2T is understood with an analogous notation for the 

The part of the action that depends only on the parti- 
cle variables Ss can be found by combining the relevant 
contributions in Eqs. (jA14|) . (|A15|) and (jA19p . 



rSsiiq,}) = TSs{{q,}Ui) ~ TSs{{q,}f=2T+l 



■zSnsiiqj] ]=T+i) 



(A23) 



Note that the exponential factors that stem from the 
bath integration measure Eq. (J18p exactly combine with 
the sums in the actions (jA16p and (|A20p . The elements 
K{j,j') of the matrix K arc easily found: 



(A19) 



Kij,j') 



1 - ki-- ■ kzT 



where So[q, in] = Sos[q] + 5os[Cri] + SosB[q, £.n] with 



1 for j ~ j' 

nUv+i ki for j > f 
for j < j' . 

(A24) 



fci ■■■ks 



R 



T-l 



T 
+ At' ^ [?i<;^.e„.j-l + gnqMh + Cnj-l)] . 

(A20) 
We introduced the imaginary time path step At' — f3h/T. 



The Gaussian integral in Eq. (jA2ip is now readily done. 
Explicit expressions of Eq. (jA24p in the continuum limit 
T — >■ 00 are easily obtained: for instance, when j, j' < T, 
K{j,j') couples to two q'^ fields and is given by 



K(j,j') = 



I _ p-fihu "l p-l3hu:+iu:At(j~j'} £qj. j <^ j' < T 

(A25) 



Note that under the sum over j and j' only its sym- 
metrized version occurs. By reintroducing the fields q^ , 
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q and q^ wc find [9| 

J-[g+,g-,(Z°]=cxp('-i$[g+,q-,5"]') , (A26) 
where the exponent reads 

- dr / Aa K{-iT + ia)q'^{T)q\a)+ dr^q° (r) 

JO Jo Jo 2 



The real and the imaginary parts of the kernel K{9) 
Kii[9) + iKi[e) are found to be 

p dc. coshK^n/2-T)] 

^ Jo n smh[ujl3h/2] ^ ' 



Kj{s — it) 



(A33) 

dw sinh[a;(/3?i/2 - r)] . 

— '5' w) . . , , sm ws) . 

TT snih[wpAi/2J 

(A34) 



/O ^0 Jo 

1/ dr / As K*{s-iT)q^{T)[q+{s)-q-{s)\ 











(A27) 



The imaginary time argument t varies in the interval 
[0,/??i] so that it is convenient to introduce the Fourier 
scries of K{s — it) with respect to r. Introducing the 
Matsubara frequencies 



At / Au [q+{s) - q-{s)] [K{s - u)q+{u) 
10 Jo 

^K*{s~u)q-{u)\ 



l^k 



27rfc 



(ASS) 



d,s^ 







1 (s) - q (s) 



we find 



The kernel K reads for complex times 6 — s — iT, < 

T<ph 



Kf 



(^-*^)=«fi E .9fc(^K'""' (A36) 



/3?i 



fc — — CO 



and 



Kio) = y: 



Qn^ cosh[a;„(/3?i/2 — i9)] 



h sinh[aj„/3?i/2] 
and the constant fi is given by 



M 



2E 



huJn 



(A28) 



(A29) 



X 



i{s-ir) = ^ ^ ^(s)e^^^^ (A37) 



/3?i 



fc— — oo 

where the functions g^ and fk are defined through 



9kis) 







°° dw „, , 2a; , , , . „„. 

S'(lj)-^ 2Cos(^s) (A38) 



TT w + i^f; 



Note that for a fcrmionic bath the only difference lies 
in the boundary conditions enforced by the trace op- 
eration: for fermions anti-periodic boundary conditions 
apply in contrast to periodic boundary conditions for 
bosons. This difference is incorporated by replacing —ki 
by fci in Eq. (|A22|) . The analysis for bosons can then be 
repeated, leading to a fermionic bath kernel where the 
cosh and the sinh in Eq. (|A28[) interchange their posi- 
tions. 

The environment can be regarded as a proper heat 
bath only if the spectrum of the harmonic oscillators be- 
comes quasi-continuous. Accordingly, we introduce the 
spectral density of the bath through 



and 



fkis) = r —S{L.)^^sm{c,s) . (A39) 

In the following we will express most quantities in terms 
of the functions gk and fk ■ For real times the real and the 
imaginary parts of the kernel (|A31[) read [see Eqs. (|A36[) 
and (1X371) ] 

Kr{s)= / — S'(cj)coth(/3nw/2)cos(ws) (A40) 

Jo T^ 

and 



2 2 i*oo J 

S{uj) = TT y ^5{uj - a;„) ^ TT y „ ''" 5{uj - (A^Q) Ki{s) = - — 5(w) sin(cjs) . (A41) 

^—'n ^-^ 2mnUJn Jo ^ 



Then the kernel K{d) and the constant fj, arc rewritten 
in terms of the spectral density 



KiO) 



Auj ^ coshM(3n/2~t9)] 
V^^""^ smh[^ph/2] ^^'^^> 



We now eliminate the local terms in Eq. (|A27[) . Wc define 
the two new kernels 

7(s) = T7/ — ^T^cosM (A42) 



and 



and 



H = 2 



Aid S{iu) 



(A32) 



AI ./o Tl" <^ 



M-) = § E a^-- 



I3h 



(A43) 



A:— — OO 
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where Cfc is defined by 



1 






1 [^ duj S{lu) 2vl 



The latter kernel is related to Kni—ir) via 

rhp 



(A44) 



/ dr / d(TA'i?(-iT + icr)/(r,cr) == (A45) 

Jo Jo 

- / dr/(r, t) + - dTdcrfc(T - cr)/(T, cr) , 



with a generic function /. In terms of the kernels 7(s), 
A:(t), A'*(s — ir), Kj^{s — u) and the linear combinations 

a; = (q+ + g")/2 and x = q+ - q^ (A46) 

the exponent of the influence functional reads 

<^[x,x,q°] = - / drdcr fc(T-cr)g°(r)qO(cr) 
2 Jo 



— i / dr dsK*{s — iT)q {t)x{s) 
Jo Jo 

1 /•* 

+ - / dsdu Knis - u)x{s)x{u) (A47) 

2 Jo 

+iA/ / d.s x(s)—- / dw 7(s — M)i(w) . 
Jo ds Jq 



Details of the derivation of Eq. (|A47P can be found in the 
thorough analysis in [9|. 



Appendix B: Classical Brownian particle in a 

harmonic potential: Initial position measurement 

and quenches in the trapping potential 

This part is meant to be a reminder on classical 
stochastic motion induced by generic baths. None of the 
results presented herein are new but they are useful to 
be confronted with the quantum results discussed in the 
body of the paper. 

The classical Brownian motion of a particle confined in 
a harmonic potential can be described by the Langevin 
equation 



q{t) + / ds j{t - s)q{s) + uj\{t) = C(i) , (Bl) 
Jo 

where ^ is a zero mean Gaussian noise |45| with corre- 
lation {£,{t)£,{s)) = -^l{\'t — s\) and with ^{t) given in 
Eq. (|42]) [l[. In the Laplace transform formulation, the 
solution to Eq. (|Bip reads 



q{\)^g+{\) |(A)+«o + (A + 7(A))g° , (B2) 



where we used Q+ defined in Eq. (|49)) and we introduced 
the initial conditions g(0) = q" and q(0) = Vq. The 
correlation function is now easily computed and it reads 

C{\,k) = {q{X)q{n)) = ±-l^±^^g+{\)g+{n) 
pM A + K 

+ g+{\)g+{n)[vl + {\ + ^{\)){K + ^{K))q'>' 

+ a+(A)g+(K) woq" [A + 7(A) + K + 7(k)] , 

(B3) 

where we used the fact that the Laplace transform of 
7(|i — s|) with respect to t and s is given by [7(A) + 
7(k)]/(A-|-k). The initial values q^ and vq can be sharp or 
drawn from a probability distribution which is typically 
of the Maxwell-Boltzmann type, that is 



/3Mowo 
Pyi ^M = -^ cxp 

ZTT 



-P 



Mo 



Mo 2 o2 
— -0? 



(B4) 



where ujq is the frequency of the initial trapping potential 
and Mo is the initial mass. From Eq. (|B4p we easily 
derive 

{q"^) = ifiMoojl)-^ and {vD ^ {(SMo)-^ . (B5) 

As long as Wo = 1^ and M = il/g the correlation function 
can be rewritten as 






ci°q(A) + ci°q(K) 

A + K 

1 7(A) + A 

(3Muj^ A2 + 7(A) + w2 



with 



(B6) 



which is the equilibrium correlation function. 

The non equilibrium correlation can be recast in the 
form 



C(A,.) = '""''\'-r'^''"''" +g.(A)g.(.) 



A + K 

/3A/a;2 Ci'=q(A)Ci<=q(K) 



2„0^ 



pMuj^q 



+ a+(A)g+(At) z;oq° [A + 7(A) + k + 7(k)] 



1 

Jm 



(B7) 



In many cases g" and vo are uncorrelated random vari- 
ables. Then Eq. ()B7[) transforms into 



1 

Jm 



C(A..), '^""W+'^;'""'> +8,(AfeW 



A + K 

-f /3A/w2ci'=q(A)Ci<=q(K) 



2/„0-^ 



l3Muj^{q 



(B8) 



Appendix C: The equilibrium initial condition 

In this Appendix we use Eq. (|54l) in the particular case 
of an equilibrium initial condition and a subsequent evo- 
lution still in equilibrium. We show how to derive the 
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equilibrium correlation function and we prove that the 
fluctuation-dissipation theorem (FDT) is satisfied with- 
out imposing time-translational invariance (TTI) as has 
been done before in the literature M- 



1. The fluctuation-dissipation theorem 

The linear response is easily found by noting that the 
external source F{s) represents a physical drift force. 
Therefore, by calculating 



m)) 



h 6 



i 5G{t) 



exp(^[^,G])|G^o 



As n{t - s)F{s) 
(CI) 



for F 7^ one finds the response function TZ{t) with re- 
spect to the external force F{t). We set e — >■ oo which 
corresponds to the absence of any initial measurement. 
By using Eq. ([M)) we obtain 



7^(t) 



1 

M 



5+{t) 



and 



^(A) 



1 



Ma;2 + A7(A)-FA2 
(C2) 



in the time and Laplace transform domains, respectively. 
These expressions are independent of the initial condi- 
tion. Therefore, the response function is the same in and 
out of equilibrium. Moreover, it is equal to the response 
function of a classical Brownian particle [4^ if it is cou- 
pled to a bath with the same friction kernel 7. 

We will confirm the validity of the FDT when the sys- 
tem is in equilibrium. We choose the initial Hamiltonian 
to be equal to the "bulk" one, that is 



w = Wo and M = Mq , 



(C3) 



so that the initial density matrix is equal to the Boltz- 
mann weight exp(— /SH) with the terms contributing to 
U given in Eqs. ([J), © and ©. The initial state is 
not perturbed by any measurement, so we take e ^ 00 
which implies A' = A and W = U. From Eq. ((57|) we 
find the equilibrium correlation function C^^(t,t') which 
in the Laplace transform version reads 

C'^'^{X,k) = A^^(A)g\(^) |aA«:+ AaCi(«;) 

+ A,C,(A) + n - ^C.(A) - ^C.(.) + ^i?"(A, .) 

(C4) 

This expression can be greatly simplified. Wc first note 
that from the definitions of gk and fk in Eqs. (|A38p 
and (|A39P it follows that 

fk{s) = ,^kgk{s) , A- (A) = ^gfc(A) , (C5) 



where we used /fc(0) ~ 0. The Laplace transform of the 
kernel 

- fk{s).fk{s')] + MoKnis - s') (C6) 

[sec Eq. ([55)1 and (^01 for A = A'] can now be written as 






dkW + gkin) 

X + K 



(C7) 



and by defining hk{X) = gk{X)/M + A we find that the 
expression in the curly brackets in the rhs of Eq. (|C4[) 
can be recast as 



Rh ^' 



^,z_.-H^-Ai)'^^(^^'^^-('^) 



-y 



pn 



hkW + hkjn) 

X + K 



(C8) 

where we used Eqs. (p6| and ([37]). By combining 
the expression for the Laplace transform of the cosine 
J^dt e"^*cos(a;f) = A/(A2 + w^) with Eqs. (1X421) 
and (jA38p wc obtain 



M 



5fc(A) 



A 



2 ,2(l'^fcl7(l^fc|)"A7(A)) 



(C9) 



By using instead Eqs. (|A42p and (jA44[) we derive 

a = I^;.i7(i^fci) ■ (CIO) 

The kernel 7 can be eliminated in favor of Q-f- through 
Eq. (gni) which yields 



hiX) 



A 



-I A2 



g^\Wk\)-g+\x) . (cii) 



This expression can now be inserted via Eq. (jCSp into 
Eq. (jC4p to find the equilibrium correlator. Note that, 
for Wo = w we have ^+(|i/fc|) = Uk- With the help of the 
one variable function 



C^"^(A) 



1 ,-^ A 

Jm^i'I-X^ 



g+{x) - g+{\iyk\) {ci2) 



the equilibrium correlation function becomes 
Ci<"i(A) -F C1<'<i(k) 



C^^iX^K) 



X + K 



(C13) 



which clearly displays time translational invariance 
(TTI). Indeed, the Laplace transform with respect to t 
and t' of a generic function f{\t — t'\) that depends only 
on the time difference is equal to [/(A) -I- /(k)]/(A -I- k), 
where /(A) denotes the Laplace transform of f{t) with 
respect to t. Hence, we haveC'"i(i,i') = C^'"i{\t-t'\)witii 
the explicit Laplace representation of C"'^'"^ in Eq. (|C12p . 
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The equilibrium correlation function is thus found with- 
out imposing TTI. By imposing TTI Eq. ((UT2)) can be 
directly found from Eq. (|57|) by setting t' — which 
simplifies the expression considerably. Remember that 
^+(i = 0) = 1 and G+{t = 0) = 0. By taking the Laplace 
transform of the result with respect to t and by using 
Eq. ([C9|) one easily recovers Eq. (|C12[) . 

It is now straightforward to establish the relation be- 
tween C^'^iit) and n{t). Firstly, we note that since C^'^'iit) 
is an even function of t its Fourier transform C^^'^{uj) is 
related to its Laplace transform through 

Ci°^(cj) = Ci°'i(iw) + Cl°'^(-^LJ) . (C14) 

Thus, by using Eqs. (IC2|) and (|C12|) we have 



C^^^itu) 






niiLu) - n{-iuj) . (C15) 



Now, since the Fourier transform of the response func- 
tion, TZ{uj), is related to its Laplace transform via TZ{uj) ~ 
TZ(iijj) due to causality we obtain the quantum FDT in 
the form 



C^°'i(a;) = hcoih[ojpn/2\lra'R{uj) 



(C16) 



where we used the formula X]fc'^/('^^ ~^ ^k) ~ 
(/3?i/2) coth[a;/3?i/2]. This result is completely general, 
in the sense that it applies to any bath, as it should. 



Appendix D: Asymptotic behavior of t/_|_(t) and 
C^'"^(t) for Ohmic dissipation 

In the case of Ohmic dissipation the spectral function 
has the form 



A ^- one finds 



Ci°^(A) 



dv A 



g+ix)-g+{iy) 



diy X^{iy^-l) + X{iyjiXiy)-j) 



v^-l 



uj-* 



(D4) 



X f°° dv i/2 _^ ^^(^) 



where the ellipsis stands for higher orders in A. Now, by 
observing the ultraviolet behavior of Eq. (|D3I) one easily 
argues that all the terms in the rhs of Eq. (ID4I) are of or- 
der ~ A. Therefore, the long time behavior of the Ohmic 
equilibrium correlation function at zero temperature is 



y-^lcq 



1 

i2 



for t ^- oo and fih >• jcj 



7' 



741" Y65) 



It is straightforward to show by direct inversion of the 
Laplace transform that the propagator Q^(t) is exponen- 
tially suppressed for large times (and for Ohmic dissipa- 
tion) on a typical time scale 7/2, hence we have 



e+(i) - e-''*/2 foj. ^ 



(D6) 



Equation (|D6|) holds for all temperatures. In the high 
temperature regime one has z/j, — >• cxd so that C^°'^(A) ~ 

- ^+(A)-l/a;2 /(/3MA). Translated into real time 

this states that ^+(t) is proportional to the derivative 
of C^'"^(i) which is nothing else than the classical FDT. 
Accordingly, we find 



^Icq 



-7*/2 



for t^oo and ^fi '^\u? - ^i^ j^Y''^) 
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SiiS) — juj for ai ^- . 



(Dl) 



For large frequencies one typically introduces a high fre- 
quency cutoff function (since the ultraviolet divergence 
is unphysical) that we choose to be of the Drude-type 
w|j/(w|) -|- a;2) where uju ^ lu is the high frequency cut- 
off. The bath kernel 
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7(i) = 7CJ_De 



-LJot 



(D2) 



then has a finite memory and a simple form in the Laplace 
domain, namely 



7(A) = 7 



ud 



ojd + X 



(D3) 



We are interested in the equilibration behavior of the 
correlation function Eq. (|62p when quantum effects dom- 
inate. In order to find the long-time behavior of C^'"^(i) 
and g+{t) we study their small-A behavior. In the very 
low temperature limit the sum over the Matsubara fre- 
quencies in Eq. (jC12p can be replaced by an integral. For 
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